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The relatively small binding energy in nuclei suggests that they may be well represented by near- 
BPS Skyrmions since their mass is roughly proportional to the baryon number A. For that purpose, 
we propose a generalization of the Skyrme model with terms up to order six in derivatives of the pion 
fields and treat the nonlinear a and Skyrme terms as small perturbations. For our special choice 
of mass term (or potential) V , we obtain well-behaved analytical BPS-type solutions with nonshell 
. configurations for the baryon density, as opposed to the more complex shell-like configurations found 

I ■ in most extensions of the Skyrme model . Along with static and (iso)rotational energies, we add to 

' the mass of the nuclei the often neglected Coulomb energy and isospin breaking term. Fitting the 

_ four model parameters, we find a remarkable agreement for the binding energy per nucleon B/A 

| > ■ with respect to experimental data. These results support the idea that nuclei could be near-BPS 

' Skyrmions. 

PACS numbers: 12.39.Dc, ll.10.Lm 
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I. INTRODUCTION 

^ n 

Q-c The idea suggested by Skyrme [l[ that baryon physics could emerge as solitons from an effective Lagrangian of 

' meson fields remains one of the most original and successful attempts for the description of the low-energy regime 
• of the theory of strong interactions (QCD). Although it predates QCD and was almost eclipsed by it, the proposal 
tIh ' gained strong support when it was realized that, in the large N c limit, QCD is equivalent to an effective theory of 
'. mesons [U, Perhaps the most important feature of the Skyrme model in that regard is that the soliton solutions 
' which arise are characterized by a conserved topological charge, the winding number, which Skyrme identified as the 
^ , baryon number. In other words, in this scheme, the baryons as well as nuclei are simply topological solitons. 

In its original formulation, the Skyrme model succeeds in predicting the properties of the nucleon within a precision 
of 30%. This is considered a rather good agreement for a two- parameter theory However, a number of general- 
izations of the model have been proposed to improve this concordance with baryon and nuclear physics. They mostly 
exploit our ignorance of the exact form of the low-energy effective Lagrangian of QCD for example, the structure of 
, the mass term [H-@], the contribution of other vector mesons @, Q or simply the addition of higher-order terms in 
derivatives of the pion fields Unfortunately, for now, QCD alone only gives hints that such extensions should 
appear and the complete determination of the effective Skyrme-like Lagrangian remains a most serious challenge. 

Despite such efforts, one of the recurring problems of Skyrme-like Lagrangians is that they almost inevitably 
lead to large binding energy for nuclei already at the classical level. A solution may be at hand by constructing 
effective Lagrangians with soliton solutions that saturate the Bogomol'nyi bound, i.e. so-called Bogomol'nyi-Prasad- 
Sommerfield type (BPS) Skyrmions, since their classical static energy grows linearly with the baryon number A (or 
atomic number) much like the nuclear mass. Support for this idea comes from a recent result from Sutcliffe |lOJ 
who found that BPS-type Skyrmions seem to emerge for the original Skyrme model when a large number of vector 
mesons are added. The additional degrees of freedom cause the mass of the soliton to decrease down to the saturation 
of the Bogomol'nyi bound. A different and more direct approach was proposed by Adam, Sanchez-Guillen, and 
Wereszczynski (ASW) [ll| by means of a prototype model consisting of only two terms: one of order six in derivatives 
of the pion fields [l2j and a second term, called the potential, which is chosen to be the customary mass term for 
pions in the Skyrme model The model leads to BPS-type compacton solutions with size and mass growing as Ai 
and A respectively, a result in general agreement with experimental observations. However, the connection between 
the ASW model and pion physics, or the Skyrme model, is more obscure due to the absence of the nonlinear a and 
so-called Skyrme terms which are of order 2 and 4 in derivatives, respectively. 

Following this picture, some of us (2l| have reexamined a more realistic generalization of the Skyrme model which 
includes terms up to order six in derivatives [l2T | in the sector where the nonlinear a and Skyrme terms are small. In 
that limit and for an appropriate choice of mass term, it is possible to find well-behaved analytical solutions for the 
static solitons. Since they saturate the Bogomol'iiyi bound, their static energy is directly proportional to A and we 
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recover some of the results in Ref. In fact, these solutions allow computing analytically the mass of the nuclei 



(static and rotational energy) in the regime where quadratic and quartic terms are small perturbations. Adjusting 
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the four parameters of the model to fit the resulting binding energies per nucleon with respect to the experimental 
data of the most abundant isotopes leads to an impressive agreement. 

These results support the idea of a BPS-type Skyrme model as the dominant contribution to an effective theory for 
the properties of nuclear matter. However, a few issues remain to be addressed before such a model is considered viable. 
We shall concentrate on two of them in this work. First, as for most extensions of the Skyrme model, the BPS-type 
models in Refs. [Tl[ and [2l[ generate shell-like configurations for the energy and baryon densities as opposed to what 
experimental data suggests, i.e. almost constant densities in the nuclei. We show here that it is possible to construct 
an effective Lagrangian which leads to nonshell configurations and still preserves the agreement with nuclear data. The 
second issue concerns the inclusion of the Coulomb energy and the isospin symmetry breaking term in the calculation 
of nuclear masses. In the context of the Skyrme model, these contributions have been thoroughly studied for A = 1 
[l^ - f20| but are usually neglected, to a first approximation, for higher A since they are not expected to overcome the 
binding energies which are usually large and also because finding the configurations is already numerically challenging 
so that only small A solutions are known (e.g. approximate toroidal, tetrahedral, cubic configurations for A = 1, 2, 3 
standard Skyrmions, respectively). However, for our type of near-BPS model, they may have a significant impact 
on the predictions given the already good agreement with data. Moreover, it turns out that the calculation of the 
Coulomb energy is simplified by the axial symmetry of the solutions and is calculable for all A. 

II. THE NEAR-BPS SKYRME MODEL 



We propose to study the model based 



with 



on the Lagrangian density 

£ = Co + £2 + £4 + £& 

-H 2 V{U) 
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where — U^d^U is the left-handed current of the meson fields represented by the SU(2) matrix U = <po + ir^i 
which obey the nonlinear condition 0q + (j>f = L 

The constants [i, a, /3, and A are left as free parameters of the model although we shall be interested in the regime 
where a and /3 are small. The original Skyrme model was built out of only the nonlinear a term, £2, and the Skyrme 
term, £4. One often adds the so-called mass term, £0, to take into account chiral symmetry breaking and generate a 
pion mass term for small fluctuations of the chiral field in V(U). Finally, the term of order six in derivatives of the 
pion fields £q is equivalent to 

£ J6 = -^B^B^ 

that was first proposed by Jackson et al. [I2j to allow for the possibility of w-meson interactions. Here, we define the 
topological (baryon) current density B^ : 

B» = --^Tr(L u L p L a ). (6) 

The boundary condition at infinity must be constant to ensure that solutions for the Skyrme field have finite energy 
but it also characterizes solutions by a conserved topological charge, 

A = J d 3 rB° = J d 3 rTr (L.LjL,) . (7) 

The static energy arising from £ 6 comes from the square of the baryon density 

^6 = f / {B\v)fd\ 
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so in a sense, it is the analog of the Coulomb energy 

E ,-\ 



4ir r - r' 



(8) 



except that instead of following the |r — r'| law, the interaction is replaced by a S— function 



Ea 



£6 
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J B° (r) B° (r') S 3 (r - r') d 3 rd 3 r'. 



In other words, the baryonic charge interacts locally. 

Historically, Cq and Cq were introduced to provide a more general effective Lagrangian than the original Skyrme 
model and indeed, the Lagrangian in (TTJ) represents the most general SU (2) model with at most two time derivatives. 
As an effective theory based on the 1/N C expansion of QCD, there no reason to believe that higher-order derivatives 
should be absent. However, since one generally relies on the standard Hamiltonian interpretation for the quantization 
procedure, higher-order time derivatives are usually avoided. 

As a result, the model has been studied extensively but remarkably, this was done only for values of parameters /x, 
a, /?, and A close to that of the original Skyrme model. Presumably these choices were made so that £ 2 and £4 would 
continue to have a significant contribution to the mass of the baryons and thereby preserve the relative successes of 
the Skyrme model in predicting nucleon properties and their link to pion physics (a is proportional to the pion decay 
constant F„). Yet this sector of the theory fails to provide an accurate description of the binding energy of heavy 
nuclei. 

Noting that this caveat may come from the fact that the solitons of the Skyrme model do not saturate the Bogo- 
mol'nyi bound, ASW proposed a model (equivalent to setting a = f) = 0) whose solutions are BPS-type solitons 
and have lower binding energies. A more realistic approach was proposed in Ref. [2l| to analyze the full Lagrangian 
([T]) in the sector where a and {3 are relatively small treating these two terms as perturbations. However, in spite of a 
very good agreement with experimental nuclear masses, there remains an obstacle to the acceptance of such model. 
Nuclear matter is believed to be uniformly distributed inside a nucleus whereas the solutions of the aforementioned 
models display shell-like configuration for the baryon and energy densities. Part of this work is to demonstrate that 
it is possible to construct an effective Lagrangian which leads to nonshell configuration and still preserves and even 
improves the agreement with nuclear mass data. 

We may write the general static solution as 



V 



in-rF 



cosF + in ■ t sinF 



where n is the unit vector 



n = (sin 9 cos $, sin 9 sin $, cos 9) 



(9) 



(10) 



Let us consider the model in ([T]) in the limit where a and /3 are small. For that purpose, we introduce the axial 
solutions for the a = f3 — case, 



F = F(r), Q = 6, $ = n 
where n is an integer. The static energy arising from [T] becomes 

9A 2 



(11) 
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J" 



dr [ fi z V + "-^—n^F' 2 



2a 
16/3 



F' 2 + (n 2 + 1) 



16 
sm 2 F 



sin 2 F 



(n 2 + 1) F' 2 + n 2 



sin 2 F 



Here F' = dF/dr and the topological charge simplifies to 

A = / F sin 2 Fdr = n 

Minimizing of the static energy for a — f3 = leads to the differential equation for F 

9X 2 n 2 sin 2 F / sin 2 F A 2 dV n 
4 2r 2 dr (— F - fl dF = °- 



(12) 



(13) 



(14) 
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The change of variable z — 2x ^^ r allows this last expression to be written in a simple form 



that can be integrated 



Regrouping the terms, we get 



BV 

sin 2 Fd z [sin 2 F (d z F)] - — = (15) 



isin 4 F(^ z ) 2 = K (16) 



sin 2 F 

dF^= =±(z- z ) (17) 
V2U 



where zq is an integration constant. Finally, the expression for F(z) can be found analytically provided the integral 
on the left-hand side is an invertible function of F. 

The potential (or so-called mass term) V in (jTTJ) is a key ingredient in the determination of the solution here. 
Unfortunately, its exact form is unknown and indeed, has been the object of several discussions 0, 0, EH- F° r 
simplicity, it is often assumed to be 

V = -Tr \U + C/ f - 2] = 1 - cos F. 

This form was considered in ASW for a = ft = in the context of BPS-Skyrmions and solving (fT7)) for F led to a 
BPS-compacton F(r) = 2 arccos (^ 1//3 r) for r £ 0, u~i , where v = j^j is a constant depending on the parameters 

A, fj,, and n. Note that f diverges as r — > and vanishes at r = 0. Since this solution saturates the Bogomol'nyi 
bound, the static energy is proportional to the baryon number A = n. 
A more general choice was introduced in Ref. 

4 

_ ^ V = Yj CkTl + C/tfc ~ 2 ] ( 18 ) 

k=l 

This form allows one to recover the chiral symmetry breaking pion mass term — im 2 7r • tt in the limit of small pion 
field fluctuations U = 2 exp(ir a 7r a /i 7 ' 7I -) and to find a relation between the pion mass and the parameter 

m 2 F 2 

J2k 2 C k = -^. (19) 
fe=i 

The case considered in Ref. [HI is a particular example of such potential with 

-C 1 =C 2 = C 3 =4C 4 = ^. (20) 
and Cfe>4 = 0. Assuming the axial solution (jlip . the potential simplifies to 

y = sin 2 ( E\ cos 6 f E\ . f2 l] 



and upon integration (|17|) . we get the solution 

F(r) = T2 arccos ( e~ vr " ) (22) 



with v = 1 ^ nX ■ Here, we use the absolute value in order to eliminate of the sign ambiguity of the arccos function. In 
order to set the baryon number to \A\ = n and the integration constant zo, we fix the boundary conditions F(0) = 
and F(oo) = =F7r for positive and negative baryon number respectively. Note that the exponential fall off of F at large 
r helps prevent some quantities such as the moments of inertia from becoming infinite. 

Unfortunately, the BPS-type models in Refs. [ll| and [2l| both lead to shell configurations for the baryon and 
energy densities which disagrees with experimental results. This is often the case for Skyrme models and it is clear 
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from expressions (|T2l and (fT31 that this behavior can be traced back to the form of the profile F(r) or more precisely 
to the derivative F'(r) which tends to zero near r = for such models. 
Let us consider the more appropriate solution of the form 

F(r) = T2 arccos (e~ ar2 ) (23) 

with a = ui and similar boundary conditions F(0) = and F(oo) = ^fir Here, since ^'(0) ^ 0, neither the baryon 
density 

„ n . . n sin F . 2an e ar / — 7 

B ( r ) = ~^^^ F =^ Vl - e- 2ar (24) 

nor the static energy density vanishes near r = 0. We find by inspection of p7|) , that this solution emerges from a 
potential similar to (|2ip . namely 

8sin 2 (f)cos^(f) 
9 In (cos 2 (f)) ' 1 J 



The logarithmic dependence in the denominator of this expression could be problematic since F = at r — and 
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F = =p7r at r = 00 but the limits for V are well defined and finite, i.e. lim r _j.o.oo V = §, respectively. It is interesting 



to note that according to (|T6|) , the square root of the potential 

3A7T / 2n sin 2 _F ,\ 3A7r 



y = — / _ — ^Lj" ) = ^lS°(r) (26) 

corresponds to the baryon radial density (the term in parenthesis) up to a multiplicative constant. Thus, in order 
to obtain a nonshell baryon density, it suffices to construct a potential V that does not vanish at small r. Such a 
potential would also imply that F'(0) 0. Our choice of potential clearly verifies this requirement but this relation 
also explains why the earlier BPS-type models could not generate a nonshell configuration, namely V ~ (1 — cos (-f-)) 
and F'(0) = in that limit. 

The expression (f25j) only applies to the axial solution (fTTj) and we need to write a more general form for V in terms 
of U if this is to be used in the expression for the Lagrangian. A simple but not unique approach to construct the 
potential is to identify cos(F/2) to the expression 

+ U + =cos 2 (jjl 
where / is the identity matrix. Then, a convenient expression for V is given by 



^ V {U) = -^Tr 
P V ' 576 



(21-U-W) (21 + u + wy 



ln((2/ + £/ + [/t)/4) 

Comparing this expression to (fT5)) allows retrieving each coefficient 

Ci = -0.129631/; 2 , C 2 = -0.100632/; 2 , C* 3 = -0.045532//, 
C 4 = -0.0105061/i 2 ,... 

such that 

J2k 2 C\ = -1.12798/i 2 = -^. 

fe=i 

Inserting expression (|23|) in (|12p , we get the static energy of the soliton in the small a and /3 approximation 



Estat — Fq + Ei + E4 + Eq 



(i 



with 

E = -\/2 (-3 + 2Vtj mr 3/2 X^ 

E 2 =U(V2- l) (n 2 + 1) + 6V2C ) tt 3 / 2 ^" 1 ^ (27) 

£ 4 = 64 (2 (l6\/2 (\/3 - l) - 15) n 2 + 2~) tt 3/2 /3i/ 1/3 

£ 6 = -a/2 (-3 + 2>/3~) ?i7r 3/2 A A i 

where ^ = sets the scale of the solution and £ is the Riemann Q— function . The terms V and Eq are proportional 
to the baryon number A = n as one expects from solutions that saturate the Bogomol'nyi bound whereas the small 
perturbations E2 = A 1 / 3 (02 + 62 A 2 ) and E4 = A" 1 / 3 (04 + b^A 2 ) have a more complex dependence. Part of this 
behavior, the overall factor A* 1 / 3 , is due to the scaling. The additional factor of A 2 comes from the axial symmetry 
of the solution Note that it is also easy to calculate analytically the root mean square radius of the baryon 

density 

rf^^'fiTS (28) 

which is consistent with experimental observation for the charge distribution of nuclei (r 2 ) 2 = r^A^ . 

In order to represent physical nuclei, we have taken into account their rotational and isorotational degrees of freedom 
and quantize the solitons. The standard procedure is to use the semiclassical quantization which is described in the 
next section. 



III. QUANTIZATION 

Skyrmions are not pointlike particles. So we resort to a semiclassical quantization method which consists in adding 
an explicit time dependence to the zero modes of the Skyrmions and applying a time-dependent (iso)rotations on the 
Skyrme fields by SU(2) matrix A(t) and B(t) 

U(r, t) = A(t)U(R(B(t))r)A(t) (29) 

where R(B(t)) is the associated SO(3) rotation matrix. The approach assumes that the Skyrmion behave as a 
rigid rotator. Upon insertion of this ansatz in the time-dependent part of the full Lagrangian ([1]), we can write the 
(iso)rotational Lagrangian as 

£rot = ->"J <•".; - aiW t jbj + ^biVijbj, (30) 

where a, = -iTtA^A and b { = iTvBB^ 

The moment of inertia tensors Uij is given by 



Ua = J d 3 r Uij = - J d 3 r [2aTr (T t Tj 



■4PTt([L pt Ti] [L pt Tj 
9A 2 
16 



^{[T u L p ] [L p , L q ] [L q , Tj] 



(31) 



where T t = iW [f ,U]. The expressions for Wij and Vij are similar except that the isorotational operator Tj is 
replaced by a rotational analog Si = —eikiXkLi as follows: 

Wv-JtrWv-JtrW^S,) (32, 
= J «/ '/• Vij = J d 3 r Uij (Tj -> Sj , T -> Si). (33) 



7 



Following the calculations in [2l[ for axial solution of the form (fTTj) , we find that all off-diagonal elements of the inertia 
tensors vanish. Furthermore, one can show that Un = U22 and £^33 can be obtained by setting n = 1 in the expression 
for Un. Similar identities hold for Vij and W%j tensors. The axial symmetry of the solution imposes the constraint 
L3 + nK 3 = which is simply the statement that a spatial rotation by an angle 9 about the axis of symmetry can be 
compensated by an isorotation of —nO about the T3 axis. It follows from expressions (f3"Tj) - (|33p that W\\ = W22 = 
for |n| > 2 and n 2 C/ 33 = nVF 33 = V33. 

The general form of the rotational Hamiltonian is given by [22T | 



L, 



22 



U22 



33 



Ek Un 

u 33 



K 2 



U22 U33 



(34) 



where (Ki) Li the body-fixed (iso)rotation momentum canonically conjugate to (ai) hi. The expression for the 
rotational energy of the nucleoli A = 1 simplifies due to the spherical symmetry 



E, 



N 



8U U 



It is also easy to calculate the rotational energies for nuclei with winding number \n\ > 2 



L K ,r,2 
1 h fif? 



(35) 



(36) 



with 



1 1 n z 

These momenta are related to the usual space-fixed isospin (I) and spin (J) by the orthogonal transformations 

h = —R(Ai)ijKj, (37) 



According to (f3"T|) and 
Hamiltonian is given by 



Ji = ~R{A2)lL j . (38) 
, we see that the Casimir invariants satisfy K 2 = I 2 and L 2 = J 2 so the rotational 



ffrot - 1 



J 2 I 2 , 

— + — +ZKl 
Vll Ull 



(39) 



We are looking for the lowest eigenvalue of H mt which depends 011 the dimension of the spin and isospin representation 
of the eigenstate \i, i^, ks)\j, j'3, ^3). For a = (3 = 0, we can show that ^ is negative and we shall assume that this 
remains true for small values of a and /3. Then, for a given spin j and isospin i, k must take the largest possible 
eigenvalue k^. Since K 2 = I 2 and L 2 = J 2 , the state with highest weight is characterized by k% = i and I3 = j and 
since nuclei are build out of A fermions we must have an isospin j < A/2. On the other hand, the axial symmetry 
of the static solutions implies that &3 = —ls/n where n = A. But for even A nuclei, k 3 must be an integer and 
\h/n\ < \j/n\ < \A/(2n)\ = 1/2 so 



< Ifcsl < 



.4 
2ri 



= 



Similarly for half-integer spin nuclei, \k$\ must be a half- integer so the only possible value is 



l<\h\< 



A 
2n 



Summarizing, if we assume for simplicity that the a and ft terms only generate small perturbations, the largest 
possible eigenvalue £3 is 



k = max(|fc3|) 



for A = even 
\ for A = odd 



(40) 
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The lowest eigenvalue of the rotational Hamiltonian H Iot for a nucleus is then given by |2l| 

+ , i(* + l) , t-2 " 



^1 



(41) 



The spin of the most abundant isotopes is fairly well known. The isospins are not so well known so we resort to 
the usual assumption that the most abundant isotopes correspond to states with lowest isorotational energy. Since 
* > [*3 1 , the lowest value that i can take is simply | «3 1 where i% = A/2 — Z. For example, the deuteron corresponds to 
A = n = 2,i = 0, j = 1, and n — 0, so the rotational energy reduces to 



E, 



D 



1 



(42) 



The explicit calculations of the rotational energy of nuclei then require only three moments of inertia which can be 
found analytically: 



U u 



T 3/2 

Is - 



24 (2V2- 1 



128 (2V2 (3 - 4\/3) (3n 2 + 1) + 3 (l2n 2 
3V2 (8^3 - 9) (3n 2 + 1) A V/ 3 



-1/3 



(43) 



V n = 



r 3/2 

W 



(2V2 



- 1 



fn 2 + 3) av- 1 



32 (32^2 (3 - 4^3) n 2 + 3 (67n 2 + 9)) /3^ 1/3 
12V2 (8V3-9) n 2 AV /3 



(44) 



and U 33 = Uu(n -> 1). 

So far, both contributions to the mass of the nucleus, -E s tat and E T0 t, are charge invariant. Since this is a symmetry 
of the strong interaction, it is reflected in the construction of the Lagrangian (p} and one expects that the two terms 
form the dominant portion of the mass. However, isotope masses differ by a few percent so this symmetry is broken 
for physical nuclei. In the next section, we consider two additional contributions to the mass, the Coulomb energy 
associated with the charge distribution inside the Skyrmion and an isospin breaking term that may be attributed to 
the up and down quark mass difference. 



IV. COULOMB ENERGY AND ISOSPIN BREAKING 



Even if we thought of a nucleus as a simple collection of individual protons and neutrons, there would be a repulsive 
electromagnetic force between protons and the process would require energy to bring these charges together. The 
result is an increase in the mass of the object by an amount corresponding to the Coulomb energy. Such an effect is 
of course also present in the Skyrmions description of nuclei since the static configuration has non- vanishing charge 
density. The electromagnetic and isospin breaking contributions to the mass have been thoroughly studied for A = 1, 
mostly in the context of the computation of the proton- neutron mass difference but are usually neglected, 

to a first approximation, for higher A since they are not expected to overcome the large binding energies predicted 
by the model. There are also practical reasons why they are seldom taken into account. The higher baryon number 
configurations of the original Skyrme model are nontrivial (toroidal shape for A = 2, tetrahedral for A = 3, etc.) 
and finding them exactly either requires heavy numerical calculations (see for example (23| ) or some kind of clever 
approximation like rational maps [24|. Moreover, the computation of the Coulomb energy is more challenging in 
general since it involves two integrations over volume. One can also argue that the Coulomb energy of Skyrmions 
is somewhat reduced by shell-like configurations of the charge densities as opposed to what it would be for a nearly 
constant spherical density found in electron scattering experiments. In our case however, we are interested in a more 
precise calculation of the nuclei masses and an estimate of the Coulomb energy is desirable, and even more so in our 
model which generates nonshell configurations. It turns out that the analytical form of the chiral angle F(r) in (f2"3"|) 
and the axial symmetry of the solution simplify the computation of the Coulomb energy. 
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Let us first consider the charge density inside Skyrmions. Following Adkins et al. we write the electromagnetic 
current 

with the baryon density and Jy 3 the vector current density, so the conserved electric charge is given by 

Z = j d\.J% M = J d\ (l& + J°?) =^ + k 3 (46) 

with ^3, the eigenvalue of third component of isospin in the body- fixed frame. The vector current is then defined as 
the sum of the left and right handed currents 

jfii _ jfj,i , jiii 
,J V ~ J R < J L 

which are invariant under SU(2)l <E> SU(2)r transformations of the form U — > LUB) . More explicitly, we get 

J^=^ J -o i -W«6 j (47) 



where Uij and Wjj are the moment of inertia densities in (|31 [) - (|33|) . In the quantized version, aj and bj are expressed 
in terms of the conjugate operators Ki and Li. Here we only need the relation 

l\, = U t3 aj - Wijbj 

Since the off-diagonal elements of Uij and Wij vanish when the solution is axially symmetric and also n 2 £/33 = nW33 = 
V33, we have 

A3 + W33&3 K3 . , 
a 3 = y = jy- + nb 3 

Inserting 03 in (|47[) . the isovector electric current density reduces to 

t03 ts U 33 

Jy = K sjr 

V33 

where U33/U33 may be interpreted here as a normalized moment of inertia density for the third component of isospin. 
Finally, the electric charge density is given by 

p(r) = >(r)+*3%^ (48) 

<5 U33 

where we have replaced k 3 by i% using the fact that the charge density of body-fixed and space-fixed frame only differs 
by a rotation. 

The Coulomb energy stored in a charge distribution p(r) takes the usual form ©. In practice, unless one considers 
very simple configurations, it is not possible to find an analytical expression for the Coulomb energy. Nonetheless, it 
is often helpful to expand p(r) in terms of normalized spherical harmonics to take care of the angular integrations 

p(r)=5>m(r)r, m '(0.0- (49) 

Following the approach described in J25|, we define the quantities 

Qlm(r) = j dr'r' l+2 p lm (r') (50) 
Jq 

which, at large distance, are equivalent to a multipole moments of the distribution. Then, each moment contributes 
to the Coulomb energy by an amount 

1 f°° 

U lm = — drr- 2l - 2 \Q lm (r)\ 2 
z ^0 Jq 
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and the total Coulomb energy associated to the distribution is given by 

oo I 



E c = E E ^ 



1=0 m=-l 



In our case, the angular dependence of the charge density is rather simple. The first part is a spherically symmetric 
contribution 



whereas the only non-trivial piece comes from the third moment density li 33 (r) and is proportional to sin 2 (8) 

U 33 = (^asin 2 F + 32/3 sin 2 F (V 2 + ™^ + ^F' 2 sin 2 F (^^)) sin 2 9 = u 33 (r) sin 2 6 
The summation (|49l) consists of only two terms 

B°(r) A^u 33 (r) 

Poo{r) = 2V7T— — + — — 1 3 

£ u 33 

4 f^ u 33 (r) . 

The expressions for moments Qoo( r ) an< i Q2o( r ) are found by integrating (|50|) analytically. Finally, we obtain the 
Coulomb energy by computing numerically the last remaining integral 



| 2 r" 4 + Ifeol 2 !- -8 ) r a * (51) 



The Coulomb energy alone cannot explain the isotope mass difference. This is particularly evident for A = 1 where 
the proton mass is known to be smaller than that of the neutron although the Coulomb energy alone would suggest 
otherwise. On the other hand, isospin is not an exact symmetry, a fact that may be traced back to the up and down 
quark mass difference. Several attempts have been made to modelize the isospin symmetry breaking term within the 
Skyrme model (l9l . [20| . Here we shall assume for simplicity that this results in a contribution proportional to the 
third component of isospin 

Ei = an 3 (52) 

with the parameter a/ fixed by setting the neutron-proton mass difference to its experimental value. Since both of 
them have the same static and rotational energies, 

AAf^ = (££ - E p c ) -ai = 1.293 MeV 

and 

ai = (E* c -E2)-AM^; 

Summarizing, the mass of a nucleus reads 

E(A,i,j,k 3 ,i 3 )=E stat (A)+E Iot (A,i,j,k 3 ) + E c (A,i 3 ) + E I (A,i 3 ) (53) 

where we have written the explicit dependence of each piece in terms of the relevant nuclear quantum numbers of the 
nuclei. The prediction depends on the parameters of the model fi, a,0, and A. 

V. RESULTS AND DISCUSSION 



The values of the parameters [i^ a, ft and A remain to be fixed. Let us first consider the case where a = j3 = 0. This 
should provide us with a good estimate for the values of /it, a, j3, and A required in the 4-parameter model ((T|) and, 
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after all, it corresponds to the limit where the minimization of the static energy leads to the exact analytical solution 



We need two input parameters to set /i and A. For simplicity, we choose the mass of the nucleon and that a nucleus 
X with zero (iso)rotational energy (i.e. a nucleus with zero spin and isospin) and neglect for now the Coulomb and 
isospin breaking energies. The total energy of these two states is according to (f2~T|) and (|4"Tj) 

E N = -V2(-3 + 2x/3V/ 2 A M + (54) 

3 V / 327r3/2V2(8V3-9)(A/x) 1/3 A 4 /3 V ' 



E x = -V2 (^-3 + 2^3) n X TT 3/2 Xid (55) 



Solving for A and fi we get 



3(-3 + 2V3) 1/4 3 3 / 4 nx 



4 (2E x ) 1/i (tt (8^3-9) (n x E N - E x )f' A 

= ((n x E N -E x )(8V3-9)f 4 ET 

M (3 7 r) 3/4 (2(2 % /3-3)) 5/4 n| 

As an example, let us examine the case where the nucleus X is Helium-4, the first doubly magic number nucleus with 
zero spin and isospin. Setting the mass of the nucleon as the average mass of the proton and neutron i.e. En = 938.919 
MeV and that of Helium-4 nucleus to E He = 3727.38 MeV, we get the numerical value A = 0.006 413 62 MeV _1 , 
a — j3 = and fj, = 14 908. MeV 2 which we shall refer as Set I. The masses of the nuclei including static, (iso)rotational, 
Coulomb, and isospin breaking contributions are then computed using (|53[) . Table I shows the relative deviation of the 
predicted with regard to experimental values of nuclear masses of a few isotopes (Set I) . The predictions are accurate 
to 0.4% or better even for heavier nuclei. Part of this accuracy is probably due to the fact that the static energy of a 
BPS-type solution is proportional to A so if it dominates, the nuclear masses should follow approximately the same 
pattern. However, the predictions remain surprisingly good for a 2-parameter model. Perhaps more relevant are the 
predictions of the binding energy per nucleon (B/A). The results are presented in Fig. [T] — Set I (solid line) and can 
be compared to the experimental values (black circles). We consider here only a subset of the table of nuclei in [2(| 
composed of the most abundant 144 isotopes. We observe a sharp rise of the binding energy per nucleon at small A 
followed by a slow linear increase for larger nuclei. The overall accuracy is of the order of 15% which is rather good 
considering the fact that the calculation involves the mass difference between the nucleus and its constituents. 
Experimentally the charge radius of the nucleus is known to behave approximately as 

with r = 1.25 fm. On the other hand, it is possible to calculate the root mean square radius for the baryon density 
[see Eq (|2"8|) ] which leads to 

{r 2 ) 1 = (2.599 fin) ^ (57) 

i 

For the charge radius (^ 2 m ) 2 , the dependence on A is more complex since it involves an additional isovector contri- 
bution (gSJ) 

where Z = 13 + A/2 is the charge of the nucleus. We get the expression 

( r v) = I drr 4 u 33 (r) 



ZU33 24/i 5 / 3 

128 (4V2 (-9 + IGVs) - 279] A 2/3 ^/i 2/3 



648 x 2 1 / 6 3 1 / 3 (72-8) ^ 4 / 3 aA 2 / 3 



13 x 2 5 / 6 3 1 / 6 (973-32) A 4 / 3 /// 3 
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where U33 also depends on A and is obtained by substituting U33 = Un(n —> 1) in (|4"3")l . Our computation verifies 
that the charge radius obeys roughly the proportionality relation 

(r c 2 m )' ~ (2.637 tm)Ai 

but overestimates the experimental value of ro by approximately a factor of 2. 

Let us now release the constraints on a and /? and allow for small perturbations from the nonlinear a and Skyrme 
term. In order to estimate the magnitude of the parameters a and /? in a real physical case, we perform two fits: Set II 
optimizes the four parameters (j,, a, /3 and A to better reproduce the masses of the nuclei while Set III tries to reach 
the best agreement with respect to the binding energy per nucleon, B/A. Both fits are performed with data from the 
same subset of the most abundant 144 isotopes as before. A summary of the results is presented in Table I while Fig. 
[1] displays the general behavior of B/A as a function of the baryon number for Sets I, II, III, and experimental values. 



Table I: Prediction versus experimental nuclear masses 





Set I 


Set II 


Set III 


Experiment 


H (10 4 MeV 2 ) 


1.490 80 


1.505 71 


1.729 55 




a (1(T 3 MeV 2 ) 





5.881 18 


22.0821 




(1(T 6 MeV ) 





-1.84877 


-5.80989 




A (10- 3 MeV" 1 ) 


6.413 62 


6.339 73 


5.536 91 




F* ( MeV) 





0.307 


0.594 


186 


m T (MeV) 




208 530 


82 300 


138 


e (10 4 ) 




-185 000 


-5380 




ro (fm) 


2.637 


2.617 


2.385 


1.23 


Nucleus X 




£ cxp (MeV) 


Nucleon 


Input 


-0.0008 


0.0020 


938.919 


2 H 


-0.0032 


-0.0048 


-0.0020 


1875.61 


3 H 


-0.0042 


-0.0057 


-0.0030 


2808.92 


4 He 


Input 


-0.0017 


-0.0009 


3727.38 


6 Li 


-0.0017 


-0.0034 


-0.0010 


5601.52 


7 Li 


-0.0014 


-0.0031 


-0.0008 


6533.83 


9 Be 


-0.0006 


-0.0023 


-0.0001 


8392.75 


io B 


-0.0004 


-0.0021 


-0.00001 


9324.44 


16Q 


0.0010 


-0.0008 


0.0009 


14 895.1 


20 Ne 


0.0010 


-0.0007 


0.0008 


18 617.7 


40 Ca 


0.0016 


0.0001 


0.0006 


37 214.7 


56p e 


0.0018 


0.0001 


0.0004 


52 089.8 


238jJ 


0.0004 


0.00001 


0.0006 


221 696 



We find that the two new sets of parameters are very close to Set I. The nonlinear a and Skyrme parameters a and 
/3 are small in magnitude but in order to make a relevant comparison, it is best to look at the relative importance of 
the contributions in ([l} and how they scale with respect to the parameters of the model, namely 





A/i 


av- 1 ^ : 




A/i 


Set I 


95.61 


: 





95.61 


Set II 


95.46 


4.408 x 10~ 5 : 


-2.255 x 10~ 4 


95.46 


Set III 


95.74 


1.418 x 10~ 4 : 


-7.904 x 10~ 4 


95.74 



for £0,^2, £4, and Cq respectively. Clearly, the nonlinear a and Skyrme terms are extremely small compared to that 
of £0 an d i- e - by at least 6 orders of magnitude. This provides support to the assumption that (f2"B"|) is a good 
approximation to the exact solution. The overall factor A/i remains approximately the same for all the sets but B/A 
turns out to be somewhat sensitive to these small variations because it involves a mass difference. Even more sensitive 
to small change in parameters is the charge radius r with 10% decrease between Set II and Set III (Table I) which 
suggests that the predicted value of ro should be taken as an estimate rather than a firm prediction. 
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FIG. 1: Binding energy per nucleon B/A as a function of the baryon number A: The experimental data (black circles) are 
shown along with predicted values (solid lines) for parametrization of Set I (a = j3 = 0) , Set II (best fit for nuclear masses), 
and Set III (best fit for B/A), respectively. 



Comparing Set II and Set III to the original Skyrme model with a pion mass term, we may identify 



and using (fT9")l we find 



m n = 1.0621-^=. 



These quantities, F w , e 2 and m„- take values (see Table I) which are orders of magnitude away for those obtained for 
the Skyrme model but this is expected since we have assumed from the start that a and /3 are relatively small. 

We find also that the Skyrme term has the wrong sign so it would destabilize the soliton against shrinking if it was 
not for the contribution of order six in derivatives which ensures stability against scale transformations. Indeed the 
term of order six was even introduced at one point to resolve some problems with this sign [l^ . In principle however, 
a negative coefficient for the Skyrme term could become problematic since the energy may no longer be bounded from 
below. One can argue that for our set of parameters, the relative weight of the £4 piece with respect to that of £0 
or Cq is so small, i.e. approximately /3i' 1 / 3 /A/i ~ — 10 -6 , and is at least partially canceled by that of the nonlinear a 
term Li so that the energy would remain bounded from below. To substantiate this point on the relative contribution 
of each term, it is useful at this point to invoke some relevant links noticed by Manton [27} between an effective SU (2) 
scalar Lagrangian and the strain tensor in the theory of elasticity. As in nonlinear elasticity theory, the energy density 
of a Skyrme field depends on the local stretching associated with the map U : R 3 H> S 3 . This is related to the strain 
tensor at a point in R 3 which is defined as 



Mij = di$dj$ where $ = (a, tt z ,tt x , tt v ) 

= _1, 

A 



r Tr[{L 4 ,L,}] 



where i,j refers to the Cartesian space coordinates. Mij is a 3 x 3 symmetric matrix with three positive eigenvalues 
A 2 , A?,, and A3. Three fundamental invariants emerges from in this simple geometrical interpretation due to 
Manton. They correspond to the Lagrangians £2, £4, and Cq and lead to the following energy densities, respectively: 

£1 = a{\\ +\\ + Xl) 

8 3 = - |/3| (\j\ 2 2 + A 2 A3 + A 2 A3) (59) 
£ 3 = 3 7 A 2 A 2 A 2 
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where we wrote for simplicity 7 = §^2-, and to the baryon density 

B° = ^^A^ (60) 

Assuming without loss of generality that Af > A| > A§, we find 

Si > SaXl £ 2 > -3 |0| XlXl £ 3 > ^XiX 2 2 Xl 

for a total energy density 

£ > tfy + £iM + 3 (\l - (a + 7 A?A*) > /i 2 F — 3 (A?A*) 

So, negative energy density contributions would come from regions where A3 < ^ ~ — 10~ 6 , in other words, where 
the baryon density B° in (foT))) is very small. Even for A§ = and a negligibly small, the energy density should be 
dominated by the potential term fi 2 V . If we consider the integrated energy density subject to the condition that 
the total baryon number A is a positive integer, then we expect the energy to be bounded from below for our set of 
parameters. 

Clearly for our axial solution, the £2 and £4 pieces of the Lagrangian do not play the same significant role in the 
stabilization of the soliton as they do in the case of the Skyrme model. The properties of the soliton are almost 
completely determined by the values of \i and A so and m, may not be so closely related to the nucleon mass scale 
as for the original Skyrme model. Perhaps the explanation for such a departure is that the parameters of the model 
are merely bare parameters and they could differ significantly from their renormalized physical values. In other words, 
we may have to consider two quite different sets of parameters: a first one, relevant to the perturbative regime for pion 
physics where F v and are close to their experimental value and a second one, that applies to the nonperturbative 
regime in the case of soliton. Unfortunately, one of the successes of the original Skyrme model is that it established a 
link between pion physics with realistic values for F w and and baryon masses. Such a link here is more obscure. 

On the other hand, the model in (TjQ) (in the regime where a and /3 are small) improves the prediction with regard 
to the properties of the nuclei of nuclear masses. Let us look more closely at the results presented in Fig. [T] These 
are in the form of the ratio of the binding energy per nucleon (B/A) as a function of the baryon (or atomic) number 
A. The experimental data (black circles) are shown along with predicted value (solid lines) for parametrization of 
Set I , Set II and Set III. Set I is the least accurate when it comes to reproducing the experimental data, especially in 
the heavy nuclei sector. Yet, the agreement remains within a 0.4% of the experimental masses which is much better 
than with the original Skyrme model. Moreover, since the ratio B/A depends on the difference between the mass of 
a nucleus and that of its constituents, it is sensitive to small variation of the nuclear masses so the results for B/A 
may be considered as rather good. The second fit (Set II), which is optimized for nuclear masses, overestimates the 
binding energies of the lightest nuclei while it reproduces almost exactly the remaining experimental values ( A > 40). 
Finally, the least square fit based on B/A (Set III) is the best fit overall but in order to provide a better representation 
for light nuclei, it abdicates some of the accuracy found in Set II for A > 40. 

This apparent dichotomy between light and heavy nuclei may be partly attributed to the (iso) rotational contribution 
to the mass. The size of nuclei grows as A$ and their moments of inertia increase accordingly. Also, the spin of the 
most abundant isotopes remains small while isospin can have relatively large values due to the growing disequilibrium 
between the number of proton and the number of neutron in heavy nuclei. Our numerical calculations reveal that the 
total effect leads to a (iso)rotational energy E Iot < 1 MeV for A > 10 for all sets of parameters considered and its 
contribution to B/A decreases rapidly as A increases. On the contrary for A < 10 the rotational energy is responsible 
for a larger part of the binding energy which means that B/A should be sensitive to the way the rotational energy 
is computed. So clearly, the shape of the baryon density will have some bearing on the predictions for the small A 
sector. 

Since part of this work is to propose a model with nonshell configuration, it is relevant to compare our result with 
a similar analysis [2l| which involves a typical shell-like configuration. For this purpose we repeated our calculation 
omitting the Coulomb and isospin breaking term. It turns out that both models are equally successful at reproducing 
data. Minimizing the square root of the mean squared deviation of B/A from its experimental value gives almost 
identical results for both models, a = 0.50 MeV per nucleon, despite generating completely different baryon and 
energy configurations. In fact, the absence of a variation in a signals somehow the equal inability for both models to 
provide an accurate description of light and heavy nuclei sectors at the same time. One could improve the agreement 
by fitting separately the parameters fi, a, (3, and A in the two sectors A > 40 and A < 40. But this would means 
introducing an arbitrary baryon number dependence on the parameter which could only be justified by introducing 
some kind of dynamical effect on /x, a, /3, and A. 
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FIG. 2: Contributions of the Coulomb and isospin effect to the binding energy per nucleon B/A as a function of the baryon 
number A for Set III: The experimental data (black circles) are shown along with contributions due to Estat + E ro t, E s t a t + 
Erot + Ec, and E to t =E sta ,t + -Erot + Eg +Ei (solid lines). 



The second motivation for this work regards this addition of the Coulomb and isospin breaking effect into the 
nuclear masses. These are often neglected in the context of the Skyrme model although they must inevitably be taken 
into account for a complete description of the nucleus. The Coulomb energy and isospin breaking term represent 
small corrections to the nuclear mass (of the order of 0.1% and 0.01%, respectively) however our results show that the 
Coulomb effect is much more significant in the calculation of the binding energy. The change in B/A is depicted in 
Fig. [2] for Set III in the separation between the red and blue lines (with and without Coulomb term, respectively). The 
effect increases almost linearly with the baryon number up to approximately 2 MeV per nucleon for the heaviest nuclei. 
It represents roughly half the Coulomb effect estimated in the Liquid Drop model (B/A) Conlomh — ~acZ(Z — l)A~ 4 ^ 3 
where the value of ac = 0.691 MeV/nucleon. On the other hand, the isospin breaking contribution due to Ej remains 
very small. Despite the magnitude of these corrections, it turns out that the optimization of the model parameters 
only yield, a slight improvement of the overall agreement with <j = 0.48 MeV per nucleon. 

To summarize, we have proposed a 4-terms model as a generalization of the Skyrme model. By choosing an 
appropriate form for the potential V, we allowed for near-BPS solitons with nonshell configurations for the baryon 
density in order to achieve a more realistic description of nuclei as opposed to the more complex configurations found 
in most extensions of the Skyrme model (e.g. A — 2 toroidal ,^4 = 3 tetrahedral, A = 4 cubic,...). Moreover, we 
introduced additional contributions to the mass of the nuclei coming from the Coulomb energy and an isospin breaking 
term. Fitting the model parameters, we find a remarkable agreement for the binding energy per nucleon B/A with 
respect to experimental data. These results suggest that nuclei could be considered as near-BPS Skyrmions. On the 
other hand, there remain some caveats. First, the Skyrme model provides a simultaneous description for perturbative 
pion interactions and nonperturbative baryons physics with single realistic values for and and baryon masses. 
The connection between the two sectors here seems to be much more intricate. Also, a much better agreement could 
be reached if one could construct a solution that would describe equally well the light and heavy nuclei. Finally, 
one would like ultimately to reproduce the observed structure of the nucleus, i.e. a roughly constant baryon density 
becoming diffuse at the nuclear surface which is characterized by a skin thickness parameter. A more appropriate 
choice of potential may be instrumental in achieving some of these goals. 

This work was supported by the National Science and Engineering Research Council of Canada. 
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